Tomographic-like reconstruction of the percolation cluster as a phase transition 
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We observed a phase transition-like behavior that is marked by the onset of the realization of the 
connectivity between two sites on a two-dimensional cross-section of a three-dimensional percolation 
cluster. This was found using contact-resistance atomic force microscopy on carbon black/polymer 
composites. The features in the current images, when presented as a function of the cut-off current, 
or as a function of the total area covered by the electrically connected objects, appear to obey a 
cluster statistics that is similar to the one predicted and observed in continuum percolation systems. 
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While the_basic theory of percolation seems to be well 
understood,!!! there are many aspects of classical perco- 
lation theory that have not been investigated intensively 
in spite of their fundamental importance and practical 
applications. In particular, there is a group of prob- 
lems associated with the fact that usually one observes 
three-dimensional (3D) structures with methods thai 
yield two-dimensional images. Such arej-jthe projection,El 
the illuminationcl and the cross-sectiona of 3D percola- 
tion systems. Conspicuous applications are astrophysicaL 
mapstl, biologicaffl-and medical scansJH meteorologicaocl 
and hydrologicatlil3 information, and the determination 
of the structure of t he . pycpl ating phase in binary in- 
homogeneous media flHEiHl3ll3 In particular, it was only 
very recently that the dynamic properties associated with 
these problems have been considered 113 This was the phe- 
nomenon of flow between two surface sites, that are con- 
nected by the underlying 3D percolation cluster. In this 
paper we report a manifestation of the above-mentioned 
problems in binary composites using an experimental 
method that depends on the dynamical aspects of flow 
within a percolation cluster. As we show below, the con- 
nection between the geometrical and dynamical aspects 
can be described by a percolation-type phase-transition- 
like behavior. 

While manifestations of|-3D properties in 2D data 
was found a long time agdl3 (on composites in which 
the electrical properties have indicated a percolation 
behaviorti't^l) the more recent development of atomic 
force microscopy (AFM) and its extensions, to stiffness 
microscopyp electrical force .microscopy (EFM)O and 
contact- resistance microscopyllj (C-AFM), have enabled 
more direct cross-section images of the percolation clus- 
ter. In particular, EFM images have been used for the. 
experimental determination of the fractal dimension,t£l 
Df, of the 3D percolating conducting phase in compos- 
ites made of carbon-black (CB) particles embedded in an 
insulating polymer 

The initial aim of our study was to reexamine the above 
EFM results using the superiorEj C-AFM method, since 
the results obtained in Ref—ll were subjected to a signif- 
icant uncertainty yieldingO a Df value in the range 2.6 
± 0.6. This permits a Df = 3 value, which is in the ho- 
mogeneous (3D) regime, and thus may shed doubt as to 



the validity of the underlying scaling assumption claimed 
in Ref. 11 (see below). The latter sijauld yield the value 
borne out from percolation theoryEH Df — 2.53. 

Our samples, with CB fillings between 10 and 23 vol%, 
were prepared from a commercial grade polyethylene and 
Vulcan XC-72 carbon blaclM This CB has 100-1000 nm- 
long aggregates that consist of fused spherical graphite 
particles that are 30r-jam in diameteiOilj. The com- 
posites were prepamdlLil in a common compounding and 
molding procedure l 19 H 2C l resulting in uniform plaques 0.25 
mm thick. The plaques were cut to long strips for the 
macroscopic measurements, described previouslyEj and 
to squares of 12x12 mm 2 for the present C-AFM study. 
A 5 mm wide silver-paint strip at the sample side served 
as a counter electrode to the local contact provided by 
the conducting tip. The C-AFM images were acquired 
more than 2 mm away from the strip, thus the current 
between the tip and the electrode is possible only if there 
is an "infinite" conducting cluster in the system, i.e., the 
current must flow between the tip and a branch of this 
infinite cluster. We call a cross section of such a branch, 
as reflected in the surface image, a "conducting island" . 

Our measurements were performed using a commer- 
cial (NT-MDT Solver) Scanning Force Microscope in the 
constant-force mode, with 0.03 N/m cantilever stiffness. 
The current maps were measured along with the topog- 
raphy for a given tip-electrode voltage, V. The silicon 
cantilever tips had curvatures of less than 35 nm, and 
were coated by a 25 nm-thick titanium-nitride film. We 
thus estimate the resolution in the topographic images to 
be better than 60 nm. The spatial resolution of the cur- 
rent maps is, however, much better, and was limited by 
our pixel size, which is of the order of 10x10 nm 2 . This 
is important considering the size of the conducting par- 
ticles (the CB aggregates, see above), since this means 
that we can assume that a pixel in an image belongs to 
a single particle. The noise level was around 0.02 nA in 
most measurements and thus in our image analysis we 
did not consider currents below 0.04 nA. The images dis- 
criminate between different current intensities (between 
and within islands) and thus enable statistical analysis 
of the current images taken with a constant voltage. In 
particular, we extract from each image the individual and 
inclusive areas of the islands for currents larger than a 



FIG. 1: 5x5 /im 2 topographic (left) and current (right) images for a sample of 15 vol% CB filling, obtained for two values of 
applied voltage, V = 1.8 V (lower scan) and V = 2.8 V (upper scan). The color scale range in the topographic images is 200 
nm, while in the current images it is 20 nA and 27 nA for the 1.8 V and 2.8 V scans, respectively. 



chosen cutoff current Iqo- 

Turning to the results, in Fig. 1 we show typical topo- 
graphic and current images obtained by the C-AFM for 
two values of V. We found no correlation between the 
surface morphology and the distribution of the " conduct- 
ing islands" . This implies that the current maps reflect 
mainly the bulk sample properties, namely, the 3D per- 
colation cluster. The increase of the number of islands 
and their individual areas with V underlines the 3D na- 
ture of the system as more percolation cluster branches 
can be detected with increasing V. The increase of the 
bias appears then to be associated with "deeper" cross 
sectioning of the conducting network of the bulk. 

To check the fractal dimension of the islands area on 
the surface, D s , we first considered samples with different 
CB contents. We measured the area covered by all the 
islands, A t , as a function of the scanned area ("the win- 
dow"), L 2 . This was done first in the regime for which 
scaling behavior is expectedjil i.e., v — 10 vol%, which 
is close to the percolation threshold, v c (= 9.3 vol%)t3, 
and L larger than the particle size (~ 0.5 fim in our 
case) but smaller than the correlation length (which, is of 
macroscopic order near the percolation thresholdJ13 say, 
larger than 10 ^m). The results of our measurements for 



various cutoff currents are shown in Fig. 2. 

For all nine cutoff currents considered (in the range 
0.04 < ico < 20 nA), the relation A t cx L Ds (expected 
from percolation theorycl) was obeyed, as seen in Fig. 2 
for four cutoff currents. The extracted exponent for each 
of the nine cutoff currents has a value within the interval 
1.35 < D s < 1.65. Our D s =X^ 0.15 exponent has 
the same value as that obtainedEltEZl from EFM measure- 
ments, but it has, as explained above, a higher certainbi 
which is very significant. Utilizing the expectatiorJiltj 
that the fractal dimension of the bulk cluster is given 
by D f = D s + 1, we find that D f = 2.50 ± 0.15, in 
agreement with the value expectedH for a 3D percolation 
cluster, 2.53. We have further established the validity of 
the above values, in addition to the smaller error bars 
and the fact that they were found for nine different cut- 
off currents, by an even more significant test. This is 
by performing the same measurements on a sample with 
v well removed from v Cl where the "normal" 3D dimen- 
sionality should be recovered. Three cutoff currents (0.5, 
4, and 20 nA) were used in a measurement of a v = 15 
vol% sample, all yielding a value of D s = 2.00 ± 0.05. 
The clear distinction between the results for the 10 vol%, 
and the 15 vol% samples is very convincing evidence 



1.000 



S 0.100 

CO 



tn 

"D 

£ 0.010 

CD 
> 

o 
o 



0.001 



Fig. 2, Shimoni (et al.) 
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FIG. 2: The measured total area covered by "conducting is- 
lands", A t , as a function of the square window size L for a 
sample with 10 vol% CB filling, for four different cutoff cur- 
rents. The lines are fits to the scaling dependence, At ~ L Ds , 
with the indicated D s values. The voltage applied in this case 
was 2V. The points in this and in other graphs represent mean 
values of data obtained for various samples and scans. 



that our surface measurements reveal properties of the 
bulk percolation cluster. This behavior appears to be 
general, as_it was found not only for the present " high- 
structure" EjE£IE3 CB/polymer composite but also for an- 
other, "low-structure", CB/polymer composite that we 
have studiedEEl, as well as for granular Ni/SiC>2 films. The 
CB system, however, is of particular interest, since it is 
not a "true" percolation system, in the sense that con- 
duction takes place via tunneling between the CB parti- 
cles even in the percolation cluster. Nevertheless, iLdoes 
conform with 3D continuum percolation behaviorO 

The results shown in Fig. 2 have significant implica- 
tions beyond the latter conclusion. We see that, system- 
atically, the lower the value of Ico the larger the value of 
A t , for the same V and L. Let us consider this signifi- 
cance. All the conducting islands on the surface belong 
to the percolation cluster, as otherwise, no current can 
flow between them and the counter-electrode. On the 
other hand, the value of Ico determines the correspond- 
ing longest possible current paths that can participate in 
the conduction (see below). We can then interpret the 
values of Ico i n terms of connectivity. The higher the Ico, 
the smaller the ensemble of conductance paths between 
the observed islands and the electrode via the percola- 
tion cluster. But then, two islands that are seen at a 
given Iqo are also connected among themselves via the 
3D percolation cluster. The lower the value of Ico, the 
longer the conductance paths included in the ensemble 
of possible paths and the larger the number of detected 
conducting islands that are interconnected via the per- 
colation cluster. Thus the variation of Ic acts as an 
electrical cross-sectioning parameter of the 3D percola- 
tion cluster. In other words, the lower the value of Ico, 
the "deeper" the (electrically) added " tomographic- like" 
section of this cluster. So, while it is obvious that all the 




Fig. 3, Shimoni et al. 

FIG. 3: (a) The dependence of the parameter C on the 
fractional total covered area, At/L 2 . Inset: the number of ob- 
served islands as a function of their size for At/L 2 = 4.6%. (b) 
The dependence of the average size of the observed "conduct- 
ing island" on the fractional total covered area, At/L 2 . The 
measurements were carried out for v — 15 vol%, L = 10/xm 
and V = 2V. The lines present fits to percolation theory-like 
expressions, see text. 



conducting islands are connected, the electrical realiza- 
tion of this connectedness is limited in our measurements 
to part of the conductance paths in the system. This in- 
terpretation gives a meaning to the Iqo ~^ limit as 
the Ico value for which all available connecting paths be- 
tween two islands are realized. We have then a transition 
from a partially realized (via the 3D percolation cluster) 
connectivity of the islands in the system to a fully realized 
connectivity of this system. Our experimental approach 
here is closely related tOptiie well-known concepts of op- 
timal and shortest paths £3 In fact, the higher the chosen 
current cutoff, the closer the collection of current paths to 
the shortest electrical (and statistically, probably also the 
geometrical "chemical distance"t3) path. Our Ico values 
then represent a hierarchy of collection of paths that de- 
viate from the shortest electrical path. The Ico — * limit 
is associated with the collection of all possible paths. 

Can the above onset of realized connectivity at Ico — * 
be described as a percolation-type phase transition? A 
priori there is a difficulty in trying to examine this ques- 
tion, since we do not know the distribution of conduc- 
tance paths-in order to apply Ico as a "thermodynamic 
parameter"E2l of a phase transition. However, there is 
an obvious parameter that Ico can be mapped to (see 
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Fig. 2) and with which we are familiar from the the- 
ory of continuum percolatiorcIlEjEJ. This parameter is 
the fractional "covered" or "occupied" area A t /L 2 that 
can be simply mapped ontO|-the site occupation proba- 
bility p in lattice percolationEJ. The choice of A t /L 2 as 
a "thermodynamic parameter" that describes the cross 
sectioning mentioned above also seems to suggest the 
value of its critical point, A tc /L 2 , as the value of A t /L 2 
when I C o — > 0. We point out again that the value of 
A tc /L 2 does not yield any particular feature in the ob- 
served image, and that the transition is "hidden" here 
by the underlying 3D percolation cluster. The parame- 
ter [A t — A tc )/L 2 may then be, as in lattice or continuum 
percolation, the parametetj-that characterizes the prox- 
imity to the critical pointElrcJ. Having this parameter, it 
was natural to determine if .cluster statistics such as in 
phase transitions iru.geperal.E3 and in continuum perco- 
lation in particular ,EjE3EJ will be observed in the present 
case. For that purpose we examined the behavior of two 
quantities, that can be readily monitored and analyzed. 
To find whether they show the behavior of a " phase tran- 
sition" we have to ch eck i f they obey the corresponding 
power-law behaviorstlE2lE3 for the same A tc - 

Following the above considerations we have plotted the 
number of observed islands, function of the is- 

land area, S, and we found that, as shown in the in- 
set of Fig. 3(a) for A t /L 2 = 4.6 %, the results can be 
well fitted by a^u-pc S~ T exp(— CS) behavior (as in per- 
colation theorjoE2l) . The specific prediction of percola- 
tion theory is that C oc (A t — A tc ) x ^ a where r and l/cr 
are dimensional-dependent exponents. As shown in Fig. 
3(a), very good fits to the data with that dependence are 
achieved for a sample of v = 15 vol%. The best fit yielded 
At c /L 2 = 13%, t = 0.25 and l/cr = 2.65. An analysis of 
different areas of this same sample yielded A tc / L 2 values 
in the range of 13 ± 1 results were obtained for other val- 
ues of v, and the corresponding values of A tc were found 
to increase with v. 

The other dependence that we evaluated was that of 
the average size of the conducting islands (the cluster size 



in percolation theoryal23) < S > on the above proximity 
parameter. < S > is defined in percolation theory as 
T,n s S 2 /Y>n s S , and is predicted to behave as (A t — A fc )~ 7 . 
We checked then whether the data obtained for < S > 
can be fitted by that power-law. The best fit for the 15 
vol% sample yielded 7 = 1.9 and, again, A tc /L 2 = 13%. 
This is shown in Fig. 3(b). 

The important conclusion that can be derived from 
Fig. 3 is that the dependencies of C and < S > on 
(A t — A t c) can be fitted very well to a power law, both 
with the same A tc and with percolation theory-like expo- 
nents. The quality of the data is not sufficient for obtain- 
ing exact values of the exponents, but the fact that the 
same A tc value is found with reasonable exponents does 
establish that we have here a behavior that can be well 
described as a phase transition. The interesting finding 
here is that while we do not have any "visible" phase 
transition in the image (such as, say, 2D connectivity), 
the quantity A t (or the hidden variable Iqo) serves as 
the "thermodynamic parameter" of the phase transition 
in the realization of the connectivity of the underlying 
3D system. The dynamic aspects, emphasized here by 
the dependence of the connectivity on V, are similar to 
the conduction under constant pressure differepee in the 
case of flow between two sites in porous media. t£l In fact, 
as will be described elsewhere, the values of A t that are 
determined for different values of V exhibit a similar be- 
havior to the one described here for different values of 
^co- 
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